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In the 'present paper general formulas for the electrical admittance of a 
piezoelectric crystal, in terms of its resonant frequencies and static and 
motional capacitances, are derived and applied to the investigation of the 
effect of electrode shape on the spectrum of resonances and. the capacitance 
ratio of the crystal. Particular attention is given in two cases of practical 
importance, namely, small piezoelectric coupling and thin crystal plates. 

I. INTRODUCTION 

Piezoelectric crystals are often used as circuit elements in filters and 
oscillators. Fig. 1 shows a typical admittance curve for such a crystal 
and Fig. 2 shows the corresponding equivalent circuit. At very low fre- 
quencies the crystal behaves like a capacitor, with a capacitance ap- 
proximately equal to the static capacitance between the driving elec- 
trodes. Due to the piezoelectric effect, an applied alternating electric 
field causes the crystal to vibrate and, at certain natural frequencies of 
free vibration, it is driven into mechanical resonance by the applied 
voltage. 

In the neighborhood of such natural frequencies, the admittance of 
the crystal is closely approximated by the simple equivalent circuit of 
Fig. A. This is the equivalent circuit commonly used in the applications. 
It is a good approximation over a frequency range proportional to the 
spacing between the resonant frequency, at which the admittance is 
infinite (in the absence of dissipation), and the antiresonant frequency, 
at which the admittance vanishes. 

For small electromechanical coupling this spacing is proportional to 
the capacitance ratio ('»/(' a . It is desirable to make this ratio as large 
as possible. Bcchmann and Parsons 1 have shown how this may be done 
in various simple cases. 
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Fig. 1 — The equivalent electrical admittance of a piezoelectric crystal. 

At high frequencies, where many resonances may occur within a nar- 
row frequency range, the simple equivalent circuit of Fig. 3 may cease 
to be applicable. However, it is always possible, at least in principle, to 
find an electrode configuration which does not excite one or more of 
these resonances. A simple example is a symmetric electrode configura- 
tion, exciting only symmetric modes of free vibration. Vormer 2 has 
shown theoretically and experimentally what the appropriate electrode 
shape is for the longitudinal vibrations of a piezoelectric bar. 

In the following, we consider both the question of capacitance ratio 
maximization and that of resonance suppression in some detail from the 
theoretical point of view. The investigation is divided into three parts. 

In Sections II through VI we consider the general problem of steady 
vibrations of a piezoelectric body. The principal tools used in these sec- 
tions are the piezoelectric analogs of various integral theorems of clas- 
sical elasticity, found, for instance, in Love. We prove, for example, 
that the material particle displacements corresponding to two different 
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Fig. 2 — The equivalent electrical circuit of a piezoelectric crystal. 
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Fig. 3 — The equivalent electrical circuit near resonance. 

modes of free vibration of a piezoelectric body are orthogonal just as 
they are for an ordinary elastic body. 

The principal results of these sections are general expressions for the 
equivalent electrical admittance of any piezoelectric body and for the 
motional capacitances. The latter expression indicates how the electrode 
configuration must be chosen to suppress a given resonance. 

Section VII is devoted to the case of small piezoelectric coupling, in 
which, to a first approximation, we may separate the mechanical and 
electrical problems completely, thus making it possible to obtain more 
explicit results than in the general case. A sample result states that the 
capacitance ratio C n /C is equal to the square of that portion of the 
input electrical energy which goes into exciting the nth mode, divided 
by the product of the electrostatic energy of the driving field and the 
strain energy of this mode of free vibration. As incidental by-products of 
our calculations, we obtain upper and lower bounds on the change in 
resonant frequency produced by a change in electrode configuration 
and an upper bound on the capacitance ratio. 

Finally, in Sections VIII and IX, we make use of the simplifications 
possible when the piezoelectric body is a thin plate, in particular obtain- 
ing an explicit relation for the electrode configuration maximizing the 
capacitance ratio. As a simple although somewhat artificial example, 
we consider the low-frequency longitudinal vibrations of a piezoelectric 
bar (the example treated by Vormer 2 ) in some detail, using this example 
to point out some of the advantages, as well as the difficulties, in the 
practical application of the foregoing techniques. 

Before commencing our discussion of these problems, a few general 
remarks may be appropriate. In order to apply the techniques presently 
proposed, considerable detailed information about the modes of free 
vibration is required. Until recently, such information was available 
only for a few special cases. However, at least for crystal plates, it ap- 
pears that such information may be obtained by the approximate theo- 
retical techniques developed by Mindlin and his co-workers* over the 

* The output of Mindlin and his coworkers in this field is so extensive that a 
complete bibliography is impossible here. Refs. 4, 5, and 6 are most closely con- 
nected with our work. 
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last decade. Experimental methods using an electric probe which may 
be moved over the surface of the crystal, such as the methods developed 
by Van Dyke and by Koga, Fukuyo, and Rhodes, 8 may also yield the 
desired information. Conversely, measurement of the electrical admit- 
tance of a given crystal for various electrode configurations and driving 
frequencies may yield considerable information about mode shapes. 

II. THE BASIC EQUATIONS 

The equations governing the steady vibrations of a piezoelectric 
crystal at angular frequency to may be written in the form 

r<y,y+ f*?Ui = 0, (1) 

Di.i = 0, (2) 

where 7\-> is the stress tensor, p the mass density, «,- the material particle 
displacement vector, and Z), the electric displacement vector. Here and 
in the following, we use Cartesian tensor notation (see, for example, 
Jeffries 9 ), in which commas denote differentiation with respect to the 
Cartesian coordinates (.i"i , .r 2 , Xa) and a repeated subscript indicates 
summation over all possible values of that subscript. Thus, for example, 
(2) is just the usual quasistatic electric field equation, stating that the 
divergence of the electric displacement vector vanishes. The symbols 
used for stress, strain, etc., are those used by Mason. 10 

In the case of a piezoelectric medium the stress 7\; and electric dis- 
placement Di are given in terms of the strain Sij and electric field E, 
by the linear anisotropic constitutive relations 

T.j = cutJStm - e kii E k , (3) 

Di = e A + eu^Stm , (4) 

where Sij and Ei are given in terms of the particle displacement w, and 
electrical potential V by the relations 

Sij = h(ui,j + «,.,), (5) 

Ei = -V,i. (6) 

The elastic constants c</*m , piezoelectric constants eitij , and dielectric 
constants e,/ satisfy the symmetry relations 

E E _ E B 

Cijkm ~ Cjikm ~ Cijmk — Ckmij > 

Gkn = e ki i , 
s s 

Hi = €ji • 
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Wc shall also have occasion to use the alternate constitutive relations 

Ta = c ijkm "S kni - h kij D k , (7) 

Ei = p ik a D k - h ikm S km , (8) 

S u = s ljkm 'T km + d kij E k , (9) 

Di = e lk r E k + d ikm T km , (10) 

where 

, , fli iovi = j, 

(0, for i * j, 

hikm = &ij Cjkm , 

D B I i 

Cijkm = Cijkm -T tlpijCpk-m , 

E E _ j j 

°ipkq Cjpmq ~ OjjO km , 

1 — > c E 

T _ .S . , 

*ij — Cij I U ipffijpq • 

The magnitude of the piezoelectric effect is specified by the piezoelectric 
coupling coefficient k, where A; 2 may be given by any one of the three 
dimensionless ratios, e 2 /ce, h 2 /cp, d'/se, with e, c, e, etc., being repre- 
sentative values of the corresponding material constants. For all real 
piezoelectric materials k is small compared with unity. 

III. BOUNDARY CONDITIONS 

The particle displacement ?/, and electric potential V, satisfying the 
preceding equations in a piezoelectric body B, are completely determined 
by the specification of certain conditions on its surface S (see Fig. 4). 
We shall always assume that the body is supported in such a fashion 
that its surface is free of tractions, i.e., 

Tum = o (ii) 




Fig. 4 — Piezoelectric body with driving electrodes. 
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on S, where w, is the outward normal vector. This assumption is not 
essential; we could consider just as well the case of a rigid (or compliant) 
support. It does, however, make the subsequent algebra somewhat 
simpler. 

The body is assumed to be driven by an alternating voltage of con- 
stant amplitude V\ , applied between electrodes So and Si , plated on 
its surface. Thus 

[Vi, on Si, 
V = (12) 

[0, on So . 

The remaining surface is assumed to be free of plating. On this por- 
tion of the surface, then, one has two conditions requiring that the po- 
tential and normal electric displacement be equal to the potential and 
normal electric displacement in the external field. We may take this 
external leakage field into account formally by simply requiring that 
the equations of the preceding section hold throughout all space (except 
on So and Si), with the stress, strain, and piezoelectric constants van- 
ishing identically outside B and the dielectric constant outside B being 
that of free space. In the following we shall indicate that an integral is 
to be taken over all space by using B x in place of B. To complete the 
set of boundary conditions, we assume that V vanishes at large distance 
from B. The external field is included here only for formal logical com- 
pleteness; in almost all practical problems it is of negligible importance. 

IV. INTEGRAL RELATIONS FOR PIEZOELECTRIC BODIES 

Suppose that u% and V are any (suitably continuous) vector and 
scalar functions, defined on B + S. Note that, by our previous con- 
vention, V is actually defined throughout all space. We shall assume 
that V is continuous across S and zero at infinity. In general, of course, 
Ui and V will not satisfy the equations of the preceding section. However, 
by introducing suitable volume and surface forces, Ft , 7\ , and volume 
and surface charges, q B , Qa , we may construct a boundary value problem 
satisfied by these functions. 

First we calculate strains and electric field components, *S,y , Ei , 
from (5) and (6) and stresses and electric displacement components 
from the constitutive relations, (3) and (4). The required distributions 
of volume force and charge are then given by the relations 

pFi = - (Tijj + pcoV), 

Qb = Di.i , 
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where co is an arbitrary positive constant. Similarly, the distributions of 
surface force and charge on S are 

1\ = Tijnj , 

qs = [Dihi , 

where n, is the outward normal to S and [Di\ is the jump in the electric 
displacement vector across <S, given by 

[D,] = (D,)ext - (Di) B . 

Now suppose we have any two pairs of such functions (it/, V), 
(u/',V") and calculate the corresponding stresses, electric displace- 
ments, etc. Then the divergence theorem yields 

f u/Tij/dB = f u/TtjnjdS - f &/T/ dB, (13) 

f V'Di," dB = - f V'[D")n i dS+ f E/D" dB, (14) 
or, in terms of the equivalent forces and charges, 

f pUi'iF/' + W "V) dB - - f u!T!' dS+ f Sh'Th" dB (15) 

Jg -Is J B 

f V'q B "dB = - f V'q s " dS + f E-D" dB. (16) 

For example, for the forced vibrations previously considered, with 
(u/,V) = {in",V") = (ui,V), co" = co, Fi" = Ti" = q B " = 0, addi- 
tion and rearrangement of (15) and (16) yields 

[ (cij^'SA + Uk a EiE k - pJuiUi) dB = ViQi , (17) 

fl OO 

where Vi is constant and the total charge on Si is 

Q x = f qs dS = f [D t ]nidS. (18) 

J Sj J s, 

Note that, by our previous convention, djkm* and w, vanish outside B. 
Equation (17) essentially states that the sum of the potential energy, 
made up of the strain energy and electrostatic energy, and the kinetic 
energy is equal to the energy passing into B through *Si due to the ap- 
plied voltage V\ . 

For free vibrations, with (u, ;, V,) = (w, n ,F",co„), (17) yields the 
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Rayleigh quotient for the natural frequency co„ , i.e., 

f (c^JSij-S^" + «fB?B?) dB 

•I J u 



f n . 

/ PUi Hi 



(11)) 



(IB 



This relation suggests that the effect of piezoelectricity is to increase the 
natural frequencies over the values they would have in an ordinary 
elastic body with the elastic constants djtm'- In Section VII we shall 
prove that this is the case for small piezoelectric coupling coefficient k. 
If we subtract (16) from (15) and rearrange terms, we obtain 

f [pu/iF/' + c/V") - V'q/] dB+f {ut'Tf - Vq/) dS 

= f (Si/Tij" - Ei'Di") dB. 

This equation still holds when we interchange primed and double-primed 
quantities. Furthermore 

Si/'Ti/ - E/'Di' = Si/Tij" - E/Di", 
and thus 



/ [pu/iF," + c/V) - v'q B "] dB+f (u!rr - V'q/) dS 

= f [ P Ui"{F! + »V) - V"q B ']dB + f (u i "T i '-V"qs')dS. 



(20) 



This is the analytical form of the so-called reciprocal theorem (see Ref. 
3, Chapter VII, p. 174) for a piezoelectric body. In words it may be 
stated as follows: 

The difference between the mechanical and electrical work done by the 
forces (including kinetic reactions) and charges of the first set, acting over 
the displacements and potential produced by the second set, is equal to the 
difference between the mechanical and electrical work done by the forces 
and charges of the second set, acting over the displacements and potential 
produced by the first. 

V. EXPANSION IN NORMAL MODES 

We now return to the problem of forced vibrations, described in Sec- 
tions I and II. The equations and boundary conditions of those sections 
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completely determine a unique particle displacement ui and potential 
V, except at certain natural frequencies w„(n = 1,2, • • • ) where a non- 
trivial solution (»/,•",]'") of the tree vibration problem (with V\ = 0) 
exists. If we denote the solution of the static problem (with a> = 0) by 
(u*,V), we may satisfy all of the preceding equations and boundary 
conditions, except (1), by setting 

m = Ui'+ f>»uA (21) 

n=\ 

V = V + E a n V", (22) 

,1 = 1 

where the a„'s are to be determined to satisfy (1). Since 

TijJ = Tij," + pw„ 2 w," = 0, 
in this case ( 1 ) is satisfied if 

00 

P ]C ( w » ~~ u)a n u" = puu'. (23) 

Now suppose we set (tt/,7') = (W\0, (ui",V") = (ui n ,V n ), with 
Wm p£ <o, t , in (20). We obtain 



w„, 2 - o>,/") f u"'u" dB = 0, (24) 



so that displacements corresponding to two different natural frequencies 
are orthogonal, just as they are in ordinary elasticity. Thus, multiplying 
(23) by Ui and integrating over B, we obtain 



P{Uk 



* — co 2 )a fc / UiufdB = pco 2 / u*Ui dB. 



Again applying (20), this time with (w/,F') = (u t k ,V k ), (mW") = 
(u°,V"), we obtain 



pco* 
Thus 



2 f Ui'ufdB = [ V'[Dt\nidS = r, f [D'budS. 

J H J S, J S l 



(I I 



co 2 r, f [D?\nidS 



(25) 



pcofr 2 (cor — co' 2 ) 
For future reference, we also note the identities 
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pco n 2 f w/V dB = / S./TV dB y (26) 

•'B ^B 

/ VADflmdS = f euiWEf - &/%') dB, (27) 

J S 1 J B 

obtained by direct application of the divergence theorem. 

VI. THE EQUIVALENT ELECTRICAL ADMITTANCE 

Using the apparatus developed in the previous sections, it is a simple 
matter to obtain a general expression for the admittance of the crystal 
body. The admittance F(a>) is the ratio of the total input current I\ 
into the crystal to the voltage Vi across its terminals. The current is the 
rate of increase of the total charge Qi on the electrode Si , where 

Qx = f [Difot dS, 

so that 

F(o>) - h/Vi = M?i/7i = juVx' 1 f [DtfoidS. 

Substituting 

[Di] = [DS] + E an[Di n ], 

n— 1 

we obtain 

F(«) = j» (c 8 + £ f C " \ , (28) 

\ n-1 CO,," — CO 2 / 

where the static capacitance C a is given by 

C. = Vf 1 f [DAm dS, (29) 

and the motional capacitances C„ by 

(f [Dflmds) 1 
C n = Vf? L . (30) 

P aj„ 2 f uSu" dB 

•'B 

If we set 

X 

Co = Ca 2-1 ^« > 
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(28) may also be written in the form 

YM = j<* (Co + E = %—) . (31) 

\ n =i 1 — or/Un / 

Fig. 2 shows the corresponding equivalent circuit with 

L„ = l/co n 2 C„, n = 1, 2, •••, 

and Fig. 1 shows the corresponding behavior of the admittance with 
frequency. Since C„ « C, for real piezoelectric materials, the admittance 
is very nearly equal to juC, , except in the vicinity of a resonant fre- 
quency, where the admittance is that of a simple series-resonant circuit, 
shunted by the capacitance C . The corresponding antiresonant fre- 
quency «„' is given approximately by 

din " _ 1 + C„ 
2 n ' 

(tin <--0 

or, since C„/Co « 1, 

CO,, CO,, 2C'o 

All of these considerations hold when the resonant frequencies are not 
too closely spaced, specifically when (co„+i — w„)/w„ is large compared 
with CJCo . At high frequencies, where many resonances may occur in 
a narrow frequency band, the region of applicability of the above simple 
single-resonance circuit may be very small, and the simple relation for 
the antiresonant frequency a>„' given by (32) may not be an adequate 
approximation. In order to circumvent this difficulty, we may choose 
the shape of the electrodes Si and S so that resonances in the vicinity 
of a desired resonance in this frequency range are not excited, i.e., so 
that the corresponding C'„'s vanish, or 



L 



[DDrndS = 0. (33) 



This is a necessary and sufficient condition that the nth resonance not 
be excited by the given plating shape. 

In order to make use of this condition, the surface charge [/),"]/i, must 
be determined, cither theoretically, perhaps by using an approximate 
theory of the type applied so successfully by Mindlin and his cowork- 
ers, 4,5 ''' or experimentally using a point electrical probe as developed by 
Van Dyke 7 or Koga, Fukuyo, and Rhodes. 8 In the case of high-frequency 
vibrations, one would expect that the resulting electrode shape would 
be quite complex. Clearly the choice of electrode shape suppressing a 
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finite number of resonances is not unique. For these and other even more 
pressing reasons, it is not clear whether the use of the condition is of any 
particular practical importance. 

Besides affecting the value of C n , the choice of electrode shape also 
affects the magnitude of the static capacitance C„ (or Co). Since C 
shunts the series-resonant elements, it is desirable to make its value as 
small as possible, i.e., to make the capacitance ratio (7,,/Co as large as 
possible at the desired resonance. Finally, because of piezoelectric cou- 
pling, the values of the resonant frequencies themselves depend on the 
electrode shape, although only to second order in the piezoelectric cou- 
pling coefficient k. In order to obtain more concrete information con- 
cerning these effects, in the next section we consider the case of small 
piezoelectric coupling. 

VII. SMALL PIEZOELECTRIC COUPLING 

In this section we shall obtain expressions for the static capacitance, 
the capacitance ratios, and the shift in natural frequencies due to a 
change in electrode shape, valid for small piezoelectric coupling coeffi- 
cient /.-. With an eye to the applications, these expressions should be the 
most convenient forms for application to the most common practical 
case, namely, a crystal plate. 

First of all, it is clear that, to first order in k, the static capacitance 
C, is simply the ordinary electrostatic capacitance in the absence of any 
piezoelectric effect. Thus, to first order, the static potential V" satisfies 
the equations 

Dt = e lk *E k s = -e,-,T/, 



in B, and 




on So . 



(In this instance, one should recall our previously agreed upon conven- 
tion concerning the treatment of the external region. ) The static capac- 
itance is then given by 



C, = Vi~* [ [DShidS = Vr 2 f tij'EiEf 

Jr. J h 



dB. (34! 



Next we consider the solution of the free vibration problem for small 
coupling. To first order in k the particle displacement u" and natural 
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frequency w n are those of the purely elastic problem, and the potential 
V" and electric displacement D" are of first order in k. Finally, the 
change in natural frequency is of second order in k. Thus we set 

Ui = Ui -f- Ui , 

V = V', 

J = co 2 (l + v"), 

where the quantities with index "0" are of zero order in k, those with 
single primes first order, and those with double primes second order in 
k. We have also temporarily dropped the index "n" in order to reduce 
the number of indices present. The governing equations and boundary 
conditions then become 

TijJ + pwoV = 0, in B, 

Ti/nj = 0, on S, 

DiJ = 0, in B x , 

V = 0, on S e = So + Si , 

Taf + m{ui" + v"u!) = 0, in B, 

Ti/'tij = 0, on S, 

where the 7\/s and D.'s remain to be defined. 

The definitions of these quantities, and thus Ui , wo, V, u", and v" 
depend on the choice of constitutive relations used. For example, if we 
use (3) and (4) for the stresses and electric displacements in terms of 
the strains and electric fields, we have 

i ij — Cijkm iJkm , 

D/ = e ik s E k ' + e ikm S km °, 

'V I' ,. F V I' „ /(' ' 

1 ij = Cijkm Okm ~ Ckij^k , 

whereas if we use (7) and (8) we have 
t ° — n D c ° 

i ij — C,jkm C>km , 

E/ = t3« s Dk' - hikmSkJ, 

T,j" = c ijkm D S km " - hkijDk'. 

Use of either of these sets of relations should give the same results (to 
the order considered) for the potential V, the electric displacement 
Di, and the natural frequency co 2 = w 2 (l + v"). We may apply the 
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reciprocal theorem, (20), to obtain a simple expression for v" in either 
case. We set {u',V') = (m,-°,0), (ui",V") = (ui",V) to obtain the 
expressions 

v B = I eifEi'E/dB/pa/ f u?u?dB, (35) 

Vd = - f PtfDi'D/ dB/pu D % f ufufdB, (36) 

J B x J B 

where u E and w D are the natural frequencies found using the elastic 
constants c,y/. m B and c ijkm D , respectively, and v E and v D are the corre- 
sponding frequency shifts due to the piezoelectric effect. Since, to second 
order, 

w 2 = leg (1 + v B ) = ud (1 + Vd) 

and v E ^ 0, v D %. 0, we must have 

oiE = (^ = us , (37) 

with equality only if v E or v D vanish. According to (35) and (36), this 
can only happen if E/ or ZV vanishes identically, i.e., if the divergence 
of the vector eikmSkm° vanishes in the former case, or if the curl of the 
vector hiitmSkm vanishes in the latter case. Thus, in general, neither of 
these bounds will be attained. We can, however, prove that the lower 
bound us is most nearly attained by a completely plated crystal and 
the upper bound a>» by a completely unplated crystal. Furthermore, 
making the electrodes larger always decreases the resonant frequencies. 
These results, which are not of central importance in our present con- 
siderations, may be obtained by the application of Schwarz's inequality 
and Green's identity to (35). 

We have now obtained expressions for the static capacitance and the 
shift in resonant frequency caused by the piezoelectric effect. We now 
derive a simple and symmetric expression for the capacitance ratio 
Cn/Cs , using (26) and (27). We have 

f VJIPfliudS = f eniWBf - Sii n E k ')dB = - f e kij S iS n E k a dB, 

J Sx J B J B 

neglecting terms of higher order in h. Also, from (26), 

pcu,,' 2 f u/'u/'dB = f S./'Tu'dB, 

which, in the present case, is either the strain energy in terms of c,-# m 
or in terms of c,-y*m' D , depending on our choice of constitutive relations. 
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To first order in h the two expressions are equal, of course. Finally, using 
(34), we obtain 






ekijEk'SifdB 



i/Ei'E/ dB)( djJ'SifSkn? (IB 



■)' 



(38) 



Direct calculation, again using Schwarz's inequality, leads to the upper 
bound 



C lt < 03 D ~ — 03b" 
(\ " 03 E' 



[39) 



Again we will not dwell on this incidental result, only remarking that it 
provides a general upper bound on a quantity which we usually wish to 
make as large as possible. An alternate form of (38), in terms of stresses, 

is 






(7 drnjEfTi? (IbJ 



»t,'B t 'E/dB\ N s^jTi/'T,,,: dBJ 



T T7T 8 771 S 



(40) 



VIII. CRYSTAL PLATES 



In this section we consider a thin piezoelectric crystal plate, driven by 
symmetrically placed electrodes on its face (.r 2 = 0, .r 2 = h), as sketched 
schematically in Fig. 5. In this case, if the thickness h is small compared 




Fig. 5 — The crystal plate. 
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with the diameter of the electrodes, the static capacitance is given ap- 
proximately by 

C s = e 22 s A/h, (41) 

where A is the electrode area. To the same approximation, the static 
Held is given by 

Ei = E, a = 0, 

[ — Vi/h, for (.ri , x%) in .4, 
E 2 = \ 

0, for (xi , xi) outside A. 

(We now neglect the external field completely.) Thus we have 



f eujESSifdB = - Vi ff etuSSixt.zJdxidxt, 



where 



8tj n (xi t x t ) = Z" 1 I 8n u (x li x li x 9 )dx» 1 
or, in terms of average stresses, 

J dkijEk'Tij" dB = - Vi fjd 2il T ij n (x 1 ,x t )dx 1 dx i , 

A 

where 

V(*ii*i) = ^ ( T ij n (x l ,x 2 ,x 3 )dx 2 . 

With 

[ et'Bi'E/dB = VfC. = e^AVf/h, 
J 11 

the capacitance ratio is given by 

S-(sz)[//* ( * i * ) * , * , I- 

A 

where 



on 



U c^JSt/'S^ dB) 



(42) 



(43) 
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in terms of strains, and 



Sijkm I ij L km dti 



P-. = 



in terms of stresses. Note that the integrals in the denominators of these 
two expressions for p» are simply normalization factors. We retain them 
only to keep the dimensions straight. 

Equation (42) gives the capacitance ratio in a form which is particu- 
larly easy to study for given p„ . We wish to know either how to choose 
A to make C n /C„ zero, in order to suppress a given resonance, or how 
to choose .4 to make C„/C» as large as possible, to increase the useful 
bandwidth of the crystal. 

Clearly the former problem has many, many solutions. If, for example, 
p„ is both positive and negative over the plate area, as will be the case 
unless we are dealing with some sort of fundamental mode, we need 
only distribute the electrode area .4 so that the integrals over the posi- 
tive and negative portions are equal. The apparent difficulty when p n 
is everywhere positive (or negative) can be circumvented by a simple 
artifice. We obtained (42) by assuming at the outset that voltages of 
the same polarity were applied between all portions of the top and bot- 
tom electrodes. If instead we imagine A to be divided into two parts, 
A + and A~, with voltages of opposite polarities applied to these parts, 
then 

// p n dxidxz = // p„dxidx 3 — II p„dxidx 3 , 

A .4+ A~ 

which clearly can be made to vanish by choosing A + and A~ so that 
the corresponding integrals are equal. 

The problem of maximizing the capacitance ratio is mostly visua- 
lized simply as a geometric problem. First of all, note that in this case we 
may replace p n by its absolute value in the integral in (42), for in re- 
gions where p„ is negative we may assume that the polarity of the driv- 
ing voltage has been reversed. This obviously will always increase the 
capacitance ratio over the value it has with unreversed polarity in such 
regions. Now we imagine the surface x% = \ pjxi , Xj) | to be erected 
above the plane .r 2 = 0. The integral over | p„ | can then be interpreted 
as the volume between this surface and the base plane, cut out by a 
cylinder of cross section A. Denoting this volume by V, we then must 
choose .4 to make V 2 /A as large as possible. This problem has a very 
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simple solution, although one which is not particularly well adapted to 
calculation. 

First of all, we observe that, for fixed area, to maximize the volume 
we should choose A to be the region bounded by the appropriate level 
curve of \p„ |. Any other choice of -4, giving the same area, must con- 
sist of this region less a subregion plus an external region having the 
same area. But the value of | p„ | for points outside the level curve is 
smaller than for points inside, so that the integral can only be decreased 
by this alternate choice. 

Let us denote the area bounded by the level curve | p„(xi ,x 3 ) | = 
p = a constant by A p and the corresponding volume by V p . We must 
now choose A „ so that V p 2 /A p is maximized, i.e., so that 

d(V;-/A p ) _ 
dA p 

Now the change in V p due to an increment A^4„ in A p is pAA p , since 
the height of the surface above the base plane is constant around its 
boundary, and thus 

d(V p 2 /A p ) _ V p (2pA p - V p ) 
dA„ " A\ 

which vanishes when pA p = \V P , i.e., when the volumes above and 
below a plane through the level curve are equal. Explicitly, to maximize 
the capacity ratio we have the condition 

pA p = = // | p n (xi ,x 3 ) | dxidx-i , (45) 

Ap 

where A p is the total area of the electrode on either plate face, bounded 
by the level curve | p n (xi , x 3 ) | = p, and p n is given by (43) or (44). 
We also assume that the electrode polarity is positive when p„ is posi- 
tive and negative when p„ is negative. Since A p must of course be less 
than or equal to the plate area, the maximum capacitance ratio given 
by (45) may not actually be attained. If this is the case, a plate wholly 
covered with electrodes gives the largest possible value to the capacitance 
ratio. It may also happen that the maximizing area A P may be partly 
bounded by a level curve and partly bounded by the plate edge. 

Bechmann and Parsons 1 have considered the excitation of piezoelec- 
tric bars and plates with partially applied electrodes and have obtained 
results similar to the above for the capacitance ratio, both theoretically 
and experimentally. 
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IX. THE LONGITUDINAL VIBRATION'S OF A PIEZOELECTRIC BAR 

To illustrate the foregoing results, we now consider a concrete example. 
To keep the analysis simple, we treat the example discussed by Vormer 2 
of the low-frequency longitudinal vibrations of the piezoelectric bar 
sketched in Fig. 6. We assume that h « a « L and restrict our discus- 
sion to vibrations of wavelength large compared with a. Then we may 
assume that all the stresses T tj are zero except T n , which we take to be 
a function of X\ only. Thus, for free vibrations, 

T u = 7V(*i), 

satisfying the equation of motion 

Tn.i" + pw„ 2 Mi" = 0, 

and the stress-strain relation 

Sn 'i E pry n 

ii — Wi,i — Sun J- n • 

Eliminating Mi" between these two equations leads to the equation 

Tu.n" + (w„/c) 2 7V = 0, 
where c = \/pSun- Together with the boundary conditions 

7V(()) = W(L) = 0, 
this equation yields 

Ti\ n (x\) = sin nirXi/L, w„ = nwc/L. 

If we assume that the width of the electrode A at x = xi/L is af(x), 
as in Fig. 6, we find that the capacitance ratio is given by 



c 



(^\[f*f(x)smnirxdx] , (46) 




L- ---J 

Fig. B — The crystal l):ir. 
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where k is the piezoelectric coupling coefficient, here given by 

Sun B t22 s ' 
A = aL is the bar-face area, and A is the electrode area, given by 



A = A f \f(x)\dx. 




We must assume that ^ \f(x) | ^ 1, and we may take into account 
portions of electrode of reversed polarity by assuming that f(x) may 
be both positive and negative. 

For the wholly covered bar (/ = 1 ), for example, (46) gives 

for n ~ 1, 3, ■ • ■ , 

for n = 2, 4, • • • . 

(Clearly an even excitation, such as the above, cannot excite the odd 
modes.) This already provides an example of resonance suppression of 
a very trivial nature. In the present case, within the limitations of the 
theory, we can actually find an electrode shape which suppresses all 
resonances except one, for the T n " ' s are orthogonal. Thus, if we set 

f(x) = sin mwx, 

we find 

(j 1 7rfc 2 /4, for n = m, 

^« (0, otherwise. 

Of course, since the application of the present simple theory is restricted 
to low frequencies, we cannot expect the above to be valid for large n 
and m. On the other hand, the low-frequency resonances are sufficiently 
widely spaced so that resonance suppression is of no particular practical 
importance in this case. Thus the present example must be taken as 
illustrative, rather than practical. 

In the present case of one-dimensional vibrations, the level lines of 
| p n \, considered in the previous section, are straight lines X\ = a con- 
stant. Thus the electrode shape maximizing the capacitance ratio con- 
sists simply of a sequence of rectangular bands of suitable widths and 
polarities, extending across the width of the bar. To determine the opti- 
mum width, it suffices to consider the fundamental mode (n = 1). For 
convenience we shift the origin of coordinates to the center of the plate 
so that 

Tn(xi) = cos (tXi/L) = cos irx, 
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and assume that the electrode band covers the interval 
must choose x according to (45), 



^ .To . We 



2a 



f x 

I cos irx dx = 2a(2x cos wxo), 



or 



tan ttTo = 2xxo , 



which is satisfied by x = 0.371, which is the value found also by Bech- 
mann and Parsons 1 in this case, i.e., 74.2 per cent of the bar length cov- 
ered by electrodes. This maximizes the capacitance ratio Ci/C, for the 
fundamental mode. The capacitance ratio for overtone modes may be 
maximized similarly, since the overtone modes may be imagined to be 
made up of a set of fundamental modes for n bars, 1/n in length, set 
end to end. 

For the fundamental mode, this maximum capacitance ratio is 0.909/if 
compared with the values 0.81 1A: 2 and 0.786A: 2 for a completely covered 
bar and for a bar with sinusoidal plating which suppresses all overtones. 
Thus we may increase the capacitance ratio about 10 per cent over its 
value for the wholly covered bar by the above technique. 

A difficulty which we do not encounter with the fundamental mode 
shows up when we consider the excitation of the first overtone (n = 2). 



(a) 



(b) 




(c) 





Fig. 7 — Various electrode configurations for excitation of first overtone in 
crystal bar. 



1280 THE BELL SYSTEM TECHNICAL JOURNAL, SEPTEMBER 1961 

In order to excite this mode, we must use electrodes having odd polarity 
about the middle of the plate, as shown in Fig. 7. For example, we may 
use complete (split) electrodes, as in Fig. 7(a), sinusoidal electrodes, as 
in Fig. 7(b), or band electrodes, as in Fig. 7(c). Especially in the first 
case, there will be a substantial contribution to the total static capaci- 
tance from the capacitance between adjoining electrodes of opposite 
polarity. This contribution has been completely neglected in our pre- 
vious calculations. To make it small compared with the ordinary plate 
capacitance, proportional to electrode area, the spacing between adjoin- 
ing electrodes of opposite polarity should be large compared with the 
plate thickness. This consideration makes the band electrodes [Fig. 7(c)], 
which maximize the capacitance ratio, particularly attractive. 

The presence of a substantial electrical field component, parallel to 
the bar faces, also changes the form of the function p n and thus the 
motional capacitance. At high frequencies, with a complicated resonance 
spectrum, this field component may excite unwanted resonances. In the 
present case of low-frequency longitudinal vibrations, this effect is 
probably not important. 
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